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Abstract. I analyse our present understanding of the Frohlich polaron with emphasis 
on the response properties, in particular optical absorption. 



I The workshop "Mott's Physics in Nanowires and Quantum Dots" provides me with an opportunity 
to recall stimulating discussions, correspondence and encounters which I had with Sir Nevill, i. a. in 
Cambridge, in Antwerpen, in Leuven, in Roma... He presented a fine scries of lectures : "Solved and 
unsolved problems for non- crystalline materials" in Antwerpen in July 1975. I enjoyed discussions 
on small polarons and hopping, on the stability of bipolarons, on electronic transport in general. 
Nevill often approached a problem deeper and deeper by a kind of Plato-type approach with ever more 
penetrating questions [H El [3] . Some aspects of the present paper are treated in more generality in the 
chapter by the present author in the book [4]. I also refer to the contributions by A. S. Alexandrov; 
S. Aubry; H. Bottger, V. V. Bryksin and T. Damker; A. M. Bratkovsky; V. Cataudella, G. De Filippis 
and C. A. Perroni; H. Fehske and S. A. Trugman; Yu. A. Firsov; M. Hohenadler and W. von der Linden; 
V. V. Kabanov; P. Kornilovitch; D. Mihailovic; A. S. Mishchenko and N. Nagaosa; Guo-meng Zhao; 
M. Zoh to Ref. [4]. See further the article by A. M. Stoneham et al. [5]. 
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1. The Frohlich polaron concept 

As is generally known, the polaron concept was introduced by Landau in 1933 [6]. Initial 
theoretical [3 [H O [TOl [11], [121 [IS] and experimental [H] works laid the fundamental 
background of polaron physics. Among the comprehensive review papers and books 
covering the subject, I refer to the reviews [15], [JG] [13, [TH] [IHl [20] . 

Significant extensions and recent developments of the polaron concept have 
been realised (see, for example, [JT] [22] [23] [Ml [25] and references therein). They 
have been invoked, e. g., to study the properties of conjugated polymers, colossal 
magnetoresistance perovskites, high-Tc superconductors, layered MgB2 superconductors, 
fuUerenes, quasi- ID conductors, semiconductor nanostructures. Polaronic phenomena 
have been clearly revealed in optical properties of semiconductor quantum wells, 
superlattices and quantum dots. 

A distinction was made between polarons in the continuum approximation where 
long-range electron-lattice interaction prevails {Frohlich, or large, polarons) [8l [TT] and 
polarons for which the short-range interaction is essential {Holstein, Holstein-Hubbard 
models). An electron or a hole, trapped by its self- induced atomic (ionic) displacement 
field in a region of linear dimension, which is of the order of the lattice constant, is 
called a small polaron [26l[27] . An excellent survey of the small-polaron physics relevant 
to the conduction phenomena in non- crystalline materials and to the metal-insulator 
transitions has been given by Mott [28] [2] . 



1.1. Scaling relations for Frohlich polarons in 2D and in 3D 

Several scaling relations connect the Frohlich polaron self-energy, the effective mass, 
the impedance Z and the polaron mobility /i in 2D to their counterpart in 3D. Those 
relations were obtained by Peeters and Devreese [15] at the level of the Feynman model, 
keeping the surface phonons for the 2D system, and are listed here: 

= Iesu (^a) , (1) 



3 V 4 

^2d(«) _ ^Ipih) (2) 

Z2D{a, fi) = ^3D (^^"' ' (3) 

where Q is the frequency of the external electromagnetic field, and 

/i2D(a) = /X3D (y") ■ 

In Eq. ([2]), ("130) and m2D ('^sd) are, respectively, the polaron- and the 

electron-band masses in 2D (3D). Expressions ([I]) to (jlj) provide a straightforward link 
between polaron characteristics in 3D with those in 2D. 

The fulfilment of the PD-scaling relation [15] is checked here for the path integral 
Monte Carlo results [46] for the polaron free energy. 
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The path integral Monte Carlo results of Ref. [IB] for the polaron free energy in 
3D and in 2D are given for a few values of temperature and for some selected values 
of a. For a check of the scaling relation, the values of the polaron free energy at 
(3 = fuvi^o/ksT = 10 are taken from Ref. [16] in 3D (Table I, for 4 values of a) 
and in 2D (Table II, for 2 values of a) and plotted in Fig. [H upper panel, with squares 
and open circles, correspondingly. (Here culo is the frequency of the long-wavelength 
longitudinal optical phonons, T is the temperature). 
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Figure 1. Upper panel: The values of the polaron free energy in 3D (squares) and 
2D (open circles) obtained by TPC'2001 gB] for /3 = 10. The data for F^d (a) are 
interpolated using a polynomial fit to the available four points (dotted line). Lower 
panel: Demonstration of the PD-scaling cf. PD'1987: the values of the polaron free 
energy in 2D obtained by TPC'2001 [46] for (3 ~ 10 (squares) are very close to the 
PD-scaled [45] values of the polaron free energy in 3D from TPC'2001 for (3—10 
(open triangles). The data for |^3Z3 (^f^) are interpolated using a polynomial fit to 
the available four points (solid line). (Reprinted with permission from [37]. ©2006, 
Societa Italiana di Fisica.) 
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In Fig. [H lower panel, the available data for the free energy from Ref [IB] are 
plotted in the following form inspired by the l.h.s. and the r.h.s parts of Eq. (QP: 
F2D (a) (squares) and f-FsD (^)(open triangles). As follows from the figure, the path 
integral Monte Carlo results for the polaron free energy in 2D and 3D very closely follow 
the PD-scaling relation of the form given by Eq. ([I]).- 



F2D (a) = 




2. Bipolarons 

When two electrons (or two holes) interact with each other simultaneously through the 
Coulomb force and via the electron-phonon-electron interaction, either two independent 
polarons can occur or a bound state of two polarons — the bipolaron — can arise 
[IHl SHI [Sni EH [52l [53] . Whether bipolarons originate or not, depends on the competition 
between the repulsive forces (direct Coulomb interaction U{r) = = |^,) and the 
attractive forces (mediated through the electron-phonon interaction). In the discussion 
of bipolarons often the ratio 

^=7^ (6) 

of the electronic and static dielectric constant is used (0 < 77 < 1). It turns out that 
bipolaron formation is favoured by smaller rj. 

The Frohlich bipolaron was analyzed [521 ES] using the Feynman path integral 
formalism. Quite analogously to the above discussed relations ([I]) to (j4]), a scaling 
relation was derived between the free energies F in two dimensions F2d{o:, U, P) and in 
three dimensions F3ci{a,U, (3): 

F2Dia,U,(3) = ^F3z.(^a,^t/,/5). (7) 

Physically this scaling relation implies that bipolaron formation will be facilitated in 2D 
as compared to 3D. 

A "phase- diagram" for the polaron — bipolaron system was introduced in Refs. [521 
153] . It is based on a generalized trial action. This phase diagram is shown in Fig. [21 For 
the 3D-case, a Frohlich coupling constant as high as 6.8 is needed to allow for bipolaron 
formation. No definite experimental evidence has been provided for the existence of 
materials with such high Frohlich coupling constant. Materials with sufficiently large 
a for Frohlich bipolaron formation in 3D might exist, but careful analysis (involving 
e. g. the study of cyclotron resonance), is needed in order to confirm this. Undoubtedly, 
some modifications to the Frohlich Hamiltonian are also necessary to describe such 
high coupling because of the more localized character of the carriers in this case which 
makes the continuum approximation less valid. The confinement of the bipolaron in 2D 
facilitates bipolaron formation at smaller a. From Fig.[2lit is seen that bipolarons in 2D 
can be stable for a > 2.9, a domain of coupling constants which is definitely realized in 
several solids. 
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Coupling constant a 

Figure 2. The phase diagram for bipolaron formation in two (dashed curve) and three 
dimensions (solid curve) is presented. Bipolarons are formed below the dashed/solid 
curves. The nonphysical part U < \f2a of the (a, J7) space is shaded. (Reprinted with 
permission from Ref. [53]. ©1991 by the American Physical Society.) 

At present no consensus exists concerning the microscopic mechanism governing the 
creation of Cooper pairs in high-Tc superconductors. Nevertheless, there is definitely 
evidence for the existence of polarons and bipolarons at least in the normal phase of 
high-Tc superconductors as manifested, e. g., in their optical absorption spectra (see 
[21] and references therein). 

3. Polaron response in 3D 

3.1. Optical absorption of a polaron at strong coupling 

The problem of the structure of the large polaron excitation spectrum constituted a 
central question in the early stages of the development of polaron theory. The exactly 
solvable polaron model of Ref. [29] was used to demonstrate the existence of the so- 
called "relaxed excited states" of large polarons [30]. In Ref. [31], and after earlier 
intuitive analysis, this problem was studied using the classical equations of motion and 
Poisson-brackets. The insight gained as a result of those investigations concerning the 
structure of the excited polaron states, was subsequently used to develop a theory of 
the optical absorption spectra of polarons. The first work was limited to the strong 
coupling limit [32]. Ref. [32] is the first work that reveals the impact of the internal 
degrees of freedom of polarons on their optical properties. 

The optical absorption of light by free Frohlich polarons was treated in Ref. [32] 
using the polaron states obtained within the adiabatic strong-coupling approximation. 
It was argued in Ref. [32], that for sufficiently large a (a > 3), the (first) relaxed 
excited state (RES) of a polaron is a relatively stable state, which can participate 
in optical absorption transitions. This idea was necessary to understand the polaron 
optical absorption spectrum in the strong-coupling regime. The following scenario of a 
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transition, which leads to a "zero-phonon" peak in the absorption by a strong-couphng 
polaron, then was suggested. If the frequency of the incoming photon is equal to 
^RES = O.OGSa^cuLO, the electron jumps from the ground state (which, at sufficiently 
large coupling, is well-characterized by "s" -symmetry for the electron) to an excited 
state ("2p"), while the lattice polarization in the final state is adapted to the "2p" 
electronic state of the polaron. In Ref. |32], considering the decay of the RES with 
emission of one real phonon, it is argued, that the "zero-phonon" peak can be described 
using the Wigner-Weisskopf formula valid when the linewidth of that peak is much 
smaller than hwi^o- 

For photon energies larger than ^res + ujj^o, a transition of the polaron towards 
the first scattering state, belonging to the RES, becomes possible. The final state of 
the optical absorption process then consists of a polaron in its lowest RES plus a free 
phonon. A "one-phonon sideband" then appears in the polaron absorption spectrum. 
This process is called one-phonon sideband absorption. The one-, two-, ... K-, ... phonon 
sidebands of the zero-phonon peak give rise to a broad structure in the absorption 
spectrum. It turns out that the first moment of the phonon sidebands corresponds to the 
Franck-Condon (FC) frequency Qpc = 0.141a;^co'Lo- To summarize, the polaron optical 
absorption spectrum at sufficiently strong coupling is characterized by the following 
features (at T = 0): 

a) An absorption peak ( "zero-phonon line" ) appears, which corresponds to a transition 
from the ground state to the first RES at ^res. 

b) For Q > Qres + ^^lo, a phonon sideband structure arises. This sideband structure 
peaks around flpc- 

Note that, at T = 0, the polaron optical absorption exhibits a zero-frequency feature 

oc 6{n). 

3.2. Optical absorption of a polaron at arbitrary coupling 

Although the optical conductivity (OC) of the Frohlich polaron model attracted 
attention for years [18], there exists no exact analytic expression for it at all coupling. 
The most successful approach is that based on the Feynman path integral technique 
as applied in [33] (DSG) and [M], where the OC is calculated starting from the 
Feynman variational model (FVM) [13] for the polaron and using Feynman's path 
integral response formalism [35] . Subsequently the path integral approach was rewritten 
in terms of the memory function formalism (MFF) [36]. These approaches are 
completely successful at small electron-phonon coupling and are able to identify key- 
excitations at intermediate and strong electron-phonon coupling. A comparison of 
the DSG results with the OC spectra derived with recently developed approximation- 
free numerical Diagrammatic Quantum Monte Carlo (DQMC) [3^ and approximate 
analytical approaches has been recently performed in Ref. |38j . 
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The polaron absorption coefficient r(0) of light with frequency VL at arbitrary 
couphng was first derived in Ref. [33] (see also [36]). It was represented in the form 



where eo is the dielectric permittivity of the vacuum and n is the refractive index of 
the medium. This expression was the starting point for a derivation of the theoretical 
optical absorption spectrum of a single Frohlich polaron at all electron-phonon coupling 
strengths in Ref. [33]. S(f2) is the so-called "memory function", which contains the 
dynamics of the polaron and depends on Q, a and temperature. The key contribution 
of the work in [33] was to introduce T{Q) in the form ([8]) and to calculate ReS(f2), 
which is essentially a (technically not trivial) Kramers-Kronig transform of the function 
ImS(f2). The function ImS(f2) had been formally derived for the Feynman polaron in 
Ref. [35] where the polaron mobility fi was found from the impedance function, i. e. the 
static limit 



The nature of the polaron excitations was analysed through this polaron optical 
absorption obtained in [331 [36]. It was confirmed in [33] that the Franck-Condon states 
for Frohlich polarons at 3 < a < 7 are nothing else but a superposition of phonon 
sidebands (of the RES). It was found in [33] that a relatively large value of the electron- 
phonon coupling strength (a > 3) is needed to stabilise the relaxed excited state of the 
polaron. At weaker coupling (0 < a < 3) only "scattering states" of the polaron, in its 
ground state, play a role in the optical absorption [33l [39] . 

In the weak coupling limit, the optical absorption spectrum ([8]) of the polaron 
is determined by the absorption of radiation energy, which is reemitted in the form 
of LO phonons. For « ^ 3, the polaron can undergo transitions toward a relatively 
stable RES (see Fig. [3]). The RES peak in the optical absorption spectrum also has a 
phonon sideband-structure, whose average transition frequency can be related to a FC- 
type transition. Furthermore, at zero temperature, the optical absorption spectrum 
of a polaron exhibits also a zero-frequency "central peak" [oc 5(f2)]. For non-zero 
temperature, this "central peak" smears out and gives rise to an "anomalous" Drude- 
type low-frequency component of the optical absorption spectrum. 

In Fig. [3] (from Ref. [33]), the main peak of the polaron optical absorption for 
a = 5 ai Q = 3.51a'Lo is interpreted as due to transitions to a RES. The "shoulder" at 
the low-frequency side of the main peak is attributed to mainly one-phonon transitions 
to polaron- "scattering states". The broad structure centered at about Q = 6.3co'lo is 
interpreted as a FC band. As seen from Fig. [3l when increasing the electron-phonon 
coupling constant to a=Q, the RES peak at Q = 4.3ciJlo stabilizes. It is in Ref. [33] 
that an all-coupling optical absorption spectrum of a Frohlich polaron, together with 
the role of RES-states, FC-states and scattering states, was ffist presented. At larger a 
(say, for a > 6), the linewidths of the absorption peaks obtained in [33] are much too 
narrow. This shortcoming was noted in [33] . 




neocmb [Q - Re^{n)f + [ImS(fi)] 



2 ' 



(8) 
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Figure 3. Optical absorption spectrum of a Frohlich polaron for a = 5 and a — Q. The 
RES peak is relatively intense compared with the FC peak at these coupling strengths. 
The frequency fi/ciJLO = u is indicated by the dashed lines. The 5- like central peaks 
(at 17 = 0) are schematically shown by vertical lines. (Reprinted with permission from 
Ref. [33]. ©1972 by the American Physical Society.) 




Recent numerical calculations of the optical conductivity for the Frohlich polaron, 
performed within a Diagrammatic Quantum Monte Carlo (DQMC) method [37], confirm 
the essential analytical results derived in Ref. [33] for a < 3. In the intermediate 
coupling regime 3 < a < 6, the low-energy behaviour and the position of the RES-peak 
in the optical conductivity spectrum of Ref. [37] follow closely the prediction of Ref. 
[33] . There are some minor differences between the two treatments in the intermediate 
coupling regime 3 < a < 6: in Ref. [37], the dominant ("RES") peak is less intense 
in the Monte-Carlo numerical simulations and the second ("FC") peak develops less 
prominently. More significant differences between the two treatments appear in the 
strong coupling regime: in Ref. [37], the dominant peak broadens and the second 
peak does not develop, giving instead rise to a flat shoulder in the optical conductivity 
spectrum at a = 6.5. These qualitative differences from the optical absorption spectrum 
of Ref. [33] can be attributed to optical processes with participation of two [IQ] or more 
phonons. The above differences can arise also due to the fact that, starting from the 
Feynman polaron model, one-phonon processes are assigned relatively more oscillator 
strength and the RES of [33] therefore tends to be more stable as compared to the 
Monte-Carlo result. 

In Fig. m DQMC optical conductivity spectrum of one polaron for a = 1 is 
compared with that obtained in Ref. [H] where a canonical-transformation formalism 
is used taking into account correlations in processes involving two LO phonons. 

The optical conductivity spectra obtained with the DQMC method [371 EHl US], 
extended memory-function formalism (MFF) [38lll3], strong-coupling expansion (SCE) 
[381 H3l m] and DSG [33], [3l] for different values of a is shown in Figs. O El The key 
results of the comparison between them are the following. First, as expected, in the 
weak coupling regime, both the extended MFF with phonon broadening [38l |13] and 




Figure 4. Onc-polaron optical conductivity Rec (51) for a = 1 calculated within the 
DQMC approach [37j (open circles) and derived using the expansion in powers of a up 
to [H] (solid curve). 

DSG [331 El] are in excellent agreement with the DQMC data [3^, showing significant 
improvement with respect to the weak coupling perturbation approach [M] which 
provides a good description of the OC spectra only for a -C 1. For 4 < a < 8, where 
DSG underestimates the peak width, the damping, introduced in the extended MFF 
approach, is significant. Comparing the peak and shoulder energies, obtained by DQMC, 
with the peak energies, given by the extended MFF, and the FC transition energies from 
the SCE, it is concluded [38] that as a increases from 6 to 10 the spectral weights rapidly 
switch from the dynamic regime, where the lattice follows the electron motion, to the 
adiabatic regime dominated by FC transitions. In the coupling range 6 < a < 10, 
adiabatic FC and nonadiabatic dynamical excitations coexist. As a > 10 the polaron 




Figure 5. Comparison of the optical conductivity calculated with the DQMC method 
(circles), the extended MFF (solid line) and DSG [321 IM) (dotted line), for four different 
values of a. (Reprinted with permission from Ref. [38|. ©2006 by the American 
Physical Society.) 
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Figure 6. a), b) and c) Comparison of tlic optical conductivity calculated with the 
DQMC method (circles), the extended MFF (sohd hue) and SCE (dashed hue) for 
three different values of a. d) The energy of the lower- and higher-frequency features 
obtained by DQMC (circles and triangles, respectively) compared (i) with the FC 
transition energy calculated from the SCE (dashed line) and (ii) with the energy of 
the peak obtained from the extended MFF (solid line). In the inset the weights of 
Franck-Condon and adiabatically connected transitions are shown as a function of a. 
We have used for ry the value 1.3. (Reprinted with permission from Ref. [S^. ©2006 
by the American Physical Society.) 



OC spectrum consists of a broad FC peak, composed of LO phonon sidebands, as 
proposed in [32], where only 1 LO phonon sideband was calculated. 

An instructive comparison between the frequency of the main peak in the optical 
absorption spectrum of continuum polarons obtained within the DSG and DQMC 
approaches for various a has been made in Ref. [H]. In Fig. [7] the energy of the 




a 



Figure 7. Main-peak positions in the DQMC optical conductivity spectra of 
continuum polarons [38j compared to the analytical DSG approach |33| . (From 
Ref. glj.) 
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main peak in the OC spectra calculated within the DSG approach [531 El] is plotted 
together with that given by DQMC [371 ESj- As seen from the figure, the main-peak 
positions, obtained within DSG, are in good agreement with the results of DQMC for 
all considered values of a, up to a = 15. The difference between the DSG and DQMC 
results is slightly larger at a = 8 and for a = 9.5, but even for those values of the 
coupling constant the agreement is close. 

In summary, the rigorous numerical results, obtained using DQMC [371 US], ^ 
modulo the linewidths for a > 6 - and the recent analytical approximations (the 
extended MFF, SCE) [38l B3l HI] confirm the essence of the mechanism for the optical 
absorption of Frohlich polarons, which we proposed in 1969 - 1971 [321 [33l 

3.3. Polaron cyclotron resonance 

The RES-like resonances in r(f2), Eq. ([8|), due to the zero's of fi— ReS(n), can effectively 
be displaced to smaller polaron coupling by applying an external magnetic field B, in 
which case the condition for a resonance becomes Vt — uOf. — ReS(f2) = {ujc = eB/nihC 
is the cyclotron frequency). Resonances in the magnetoabsorption governed by this 
condition have been clearly observed and analysed. 

The results for the polaron optical absorption [33l [34] paved the way for an all- 
coupling path-integral based theory of the magneto-optical absorption of polarons (see 
Ref. [55]). This work was motivated i. a. by the insight that magnetic fields can 
stabilise relaxed excited polaron states, so that information on the nature of relaxed 
excited states might be gained from the cyclotron resonance of polarons. 

Some of the subsequent developments in the field of polaron cyclotron resonance 
are discussed below. 

Evidence for the polaron character of charge carriers in AgBr and AgCl was obtained 
through high-precision cyclotron resonance experiments in magnetic fields up to 16 T 
(see [56]). The quantitative interpretation of the cyclotron resonance experiment in 
AgBr and AgCl [56] by the theory of Ref. [55] provided one of the most convincing and 
clearest demonstrations of Frohlich polaron features in solids. 

The analysis in Ref. [56] leads to the following polaron coupling constants: a = 1.53 
for AgBr and a = 1.84 for AgCl. The corresponding polaron masses are: m* = 0.2818me 
for AgBr and m* = 0.3988me for AgCl {rrig is the mass of the free electron). For most 
materials with relatively large Frohlich coupling constant the band mass is not 
known. The study in Ref. [56] is an example of the detailed analysis of the cyclotron 
resonance data that is necessary to obtain accurate polaron data like a and rrib for a 
given material. 

Early infrared-transmission studies of hydrogen-like shallow-donor-impurity states 
in n-CdTe were reported in Ref. [57]. By studying the Zeeman splitting of the 
(Is — * 2p,m = ±1) transition in the Faraday configuration at magnetic fields up to 
~ 16 T, a quantitative determination of polaron shifts of the energy levels of a bound 
electron was performed. The experimental data were shown to be in fair agreement with 
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the weak-coupling theory of the polaron Zeeman effect. In this comparison, however, 
the value a = 0.4 had to be used instead of a = 0.286, which would follow from the 
definition of the Frohlich coupling constant. Similarly, the value a ~ 0.4 was suggested 
(see Ref. [58]) for the explanation of the measured variation of the cyclotron mass with 
magnetic field in CdTe. This discrepancy gave rise to some discussion in the literature 
(see, e. g., Refs. [591 [60] and references therein). 

In Ref. [61], far-infrared photoconductivity techniques were applied to study the 
energy spectrum of shallow In donors in CdTe layers and experimental data were 
obtained on the magnetopolaron effect. An overall good agreement was found between 
experiment and a theoretical approach, where the electron-phonon interaction is treated 
within a second-order improved Wigner-Brillouin perturbation theory and a variational 
calculation is performed for the lowest-lying donor states (Is, 2p^, 2s, 2p^, Sc?"^, 4/~^). 
This agreement was obtained with the coupling constant a = 0.286. 

The energy spectra of polaronic systems such as shallow donors ("bound polarons"), 
e. g., the Do and D~ centres, constitute the most complete and detailed polaron 
spectroscopy realised in the literature, see for example Fig. [81 

4. Many-polaron systems in 3D and in 2D 
4-1- Optical absorption of many-polaron systems 

In Ref. [6l], starting from the many-polaron canonical transformations and the 
variational many-polaron wave function introduced in Ref. [65], the optical absorption 
coefficient of a many-polaron gas has been derived. The real part of the optical 
conductivity of the many-polaron system was obtained in an intuitively appealing form 



where no is the density of charge carriers, Vk is the electron-phonon interaction amplitude 
and is the z-component of the wave vector. Equation Qj is rigorous to order 
a. This approach to the many-polaron optical absorption allows one to include the 
many-body effects in terms of the dynamical structure factor S(k, Q — culo) of the 
electron (or hole) system (e. g. in the RPA-approximation) . The experimental peaks 
in the mid-infrared optical absorption spectra of cuprates (Fig. [9]) and manganites 
(Fig. [TOl) can be adequately interpreted within this theory. As seen from Fig.[Tni for the 
case of La2/3Sri/3Mn03 the many-polaron approach describes the experimental optical 
conductivity better than the methods of [MII67]. In Ref. [68], the interesting possibility 
of coexistence of small and large polarons in the same solid was suggested. 

4-2. Ripplopolarons in multi- electron bubbles in liquid helium 

An interesting 2D system consists of electrons on films of liquid He [701 [71] . In this 
system the electrons couple to the ripplons of the liquid He, forming "ripplopolarons". 




k 



(9) 
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Figure 8. The Is 2p^,2pz transition energies as a function of magnetic field 
for a donor in GaAs. The theoretical results are compared to the experimental data 
(solid dots) [63] for the following cases: (a) without the effect of polaron and band 
nonparabolicity (thin dashed curves); (b) with polaron correction (dotted curves); (c) 
including the effects of polaron and band nonparabolicity (solid curves). (Reprinted 
with permission from Ref. [62j . (c)1993 by the American Physical Society.) 

The effective coupling can be relatively large and self-trapping can result. The acoustic 
nature of the ripplon dispersion at long wavelengths induces the self-trapping. 

Spherical shells of charged particles appear in a variety of physical systems, such 
as fuUerenes, metallic nanoshells, charged droplets and neutron stars. A particularly 
interesting physical realization of the spherical electron gas is found in multielectron 
bubbles (MEBs) in liquid helium-4. These MEBs are 0.1 fim - 100 fim sized cavities 
inside liquid helium, that contain helium vapour at vapour pressure and a nanometer- 
thick electron layer, anchored to the surface of the bubble [721, [73]. They exist as a result 
of equilibrium between the surface tension of liquid helium and the Coulomb repulsion 
of the electrons [TU [75]. Recently proposed experimental schemes to stabilize MEBs 
[76] have stimulated theoretical investigation of their properties (see e. g. [77]). 

The dynamical modes of an MEB were described by considering the motion of the 
helium surface ("ripplons") and the vibrational modes of the electrons together. In 
particular, the case when the ripplopolarons form a Wigner lattice was analyzed in Ref. 
[78] . The interaction energy between the ripplons and the electrons in the multielectron 
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bubble is derived from the following considerations: (i) the distance between the layer 
electrons and the helium surface is fixed (the electrons find themselves confined to 
an effectively 2D surface anchored to the helium surface) and (ii) the electrons are 
subjected to a force field, arising from the electric field of the other electrons. To 
study the ripplopolaron Wigner lattice at nonzero temperature and for any value of the 
electron-ripplon coupling, the variational path-integral approach [13] is used. 

In their treatment of the electron Wigner lattice embedded in a polarizable Medium, 
such as a semiconductor or an ionic solid, Fratini and Quemerais [79] described the effect 
of the electrons on a particular electron through a mean-field lattice potential. The 
(classical) lattice potential Viat is obtained by approximating all the electrons acting on 
one particular electron by a homogenous charge density in which a hole is punched out; 
this hole is centered in the lattice point of the particular electron under investigation 
and has a radius given by the lattice distance d. 

The Lindemann melting criterion [80l [8T| states in general that a crystal lattice of 
objects (be it atoms, molecules, electrons, or ripplopolarons) will melt when the average 
displacement of the objects from their lattice site is larger than a critical fraction 5q of 
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Figure 9. The infrared optical absorption of Nd2Cu04_5 {5 < 0.004) as a function of 
frequency. The experimental results of Ref. [66] are presented by the thin full curve. 
The experimental 'd-band' is clearly identified, rising in intensity at about 600 cm~^, 
peaking around 1000 cm~^, and then decreasing in intensity above that frequency. 
The dotted curve shows the single polaron result calculated according to Ref. [33] . 
The bold full curve presents the theoretical results of Ref. [64] for the interacting 
many-polaron gas with ng = 1.5 x 10^^ cm~^. a = 2.1 and m = O.Stoq. (Reprinted 
with permission from Ref. [64]. ©2001 by the American Physical Society.) 
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Figure 10. Comparison of the measured mid-infrared optical conductivity in 
La2/3Sri/3Mn03 at T = 6 K to that given by several model calculations for 
m = Swe, a of the order of 1 and rig 6 x 10^^ cm~'^. The onc- 

polaron limit [the weak-coupling approach by V. L. Gurcvich. I. G. Lang, and 
Yu. A. Firsov [54| (GLF model) and the phenomcnological approach by D. Emin 
[67j (E model)] lead to narrower polaron peaks than a peak with maximum at 
~ 900 cm~^ given by the many-polaron treatment by J. Tempere and J. T. Devreese 
(TD model) of Ref. [M]- (Reprinted with permission after Ref. [69]. ©2004 by the 
American Physical Society.) 



the lattice parameter d. It would be a strenuous task to calculate, from first principles, 
the exact value of the critical fraction but for the particular case of electrons on a 
helium surface, we can make use of an experimental determination. Grimes and Adams 
[82] found that the Wigner lattice melts when F = 137 ± 15, where F is the ratio 
of potential energy to the kinetic energy per electron. At temperature T the average 
kinetic energy of an electron in a lattice potential Viat-, characterized by the frequency 
parameter ujiat-, is 

Eun = ^^Ot\.(^\, (10) 



2kBT^ 

and the average distance that an electron moves out of the lattice site is determined by 



(r^) = coth ^ = (11) 

rrieUJiat \2kBT J rrteWf;^ 

From this one finds that for the melting transition in Grimes and Adams' experiment 
[82] . the critical fraction equals 6q ~ 0.13. This estimate is in agreement with previous 
(empirical) estimates yielding 5o ~ 0.1 [83]. 

Within the approach of Fratini and Quemerais [79] , the Wigner lattice of (ripplo)- 
polarons melts when at least one of the two following Lindemann criteria are met: 

6r = y^^3. > (12) 
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dp = — - — > do- (13) 

where p and Rcms are, respectively, the relative coordinate and the center of mass 

coordinate of the model system: if r is the electron coordinate and R is the position 

coordinate of the fictitious ripplon mass M, they are 

^ Tiler + MR ^ , 

^cms = ; P = r - R. 14 

rrie + M 

The appearance of two Lindemann criteria takes into account the composite nature of 
(ripplo)polarons. As follows from the physical meaning of the coordinates p and Rcms, 
the first criterion (fT2|) is related to the melting of the ripplopolaron Wigner lattice 
towards a ripplopolaron liquid, where the ripplopolarons move as a whole, the electron 
together with its dimple. The second criterion f|T3|) is related to the dissociation of 
ripplopolarons: the electrons shed their dimple. 

The path-integral variational formalism allows us to calculate the expectation values 
(R^^^) and (p^) with respect to the ground state of the variationally optimal model 
system. 

Numerical calculation shows that for ripplopolarons in an MEB, the inequality 

(Kms) « {P') (15) 

is realized. As a consequence, the destruction of the ripplopolaron Wigner lattice in a 
MEB occurs through the dissociation of ripplopolarons, since the second criterion (fT3l) 
will be fulfilled before the first ( fT2l) . The results for the melting of the ripplopolaron 
Wigner lattice are summarized in the phase diagram shown in Fig. [TTl 

For any value of A^, pressure p and temperature T in an experimentally accessible 
range. Fig. [H] shows whether the ripplopolaron Wigner lattice is present (points above 
the surface) or the electron liquid (points below the surface). Below a critical pressure 
(on the order of 10"^ Pa) the ripplopolaron solid will melt into an electron liquid. This 
critical pressure is nearly independent of the number of electrons (except for the smallest 
bubbles) and is weakly temperature dependent, up to the hehum critical temperature 
5.2 K. This can be understood since the typical lattice potential well in which the 
ripplopolaron resides has frequencies of the order of THz or larger, which correspond to 
~ 10 K. 

The new phase that was predicted in [78], the ripplopolaron Wigner lattice, will 
not be present for electrons on a flat helium surface. At the values of the pressing field 
necessary to obtain a strong enough electron-ripplon coupling, the fiat helium surface 
is no longer stable against long- wavelength deformations [H]. Multielectron bubbles, 
with their different ripplon dispersion and the presence of stabilizing factors such as 
the energy barrier against fissioning [85], allow for much larger electric fields pressing 
the electrons against the helium surface. The regime of A^, p, T parameters suitable 
for the creation of a ripplopolaron Wigner lattice lies within the regime that would be 
achievable in recently proposed experiments, aimed at stabilizing multielectron bubbles 
[76] . The ripplopolaron Wigner lattice and its melting transition might be detected by 
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Figure 11. The phase diagram for the spherical 2D layer of electrons in a MEB. Above 
a critical pressure, a ripplopolaron solid (a Wigner lattice of electrons with dimples 
in the helium surface underneath them) is formed. Below the critical pressure, the 
ripplopolaron solid melts into an electron liquid through dissociation of ripplopolarons. 
(Reprinted with permission from Ref. [T^- ©2003, EDP Sciences, Societa Italiana di 
Fisica, Springer.) 



spectroscopic techniques j82| [86] probing for example the transverse phonon modes of 
the lattice [87]. 

5. Polarons in 2D and in quasi-2D structures 

5.1. Polarons and cyclotron resonance in quantum wells 

The resonant magnetopolaron coupling [88] in bulk manifests itself near the LO-phonon 
frequency for low electron densities (see, e. g., Refs. [89], [901 EI]) and also for higher 
electron densities |92j . 

For example, Nicholas et al. [H9] demonstrated polaron coupling effects using 
cyclotron resonance measurements in a 2DEG, which naturally occurs in the polar 
semiconductor InSe. One clearly sees, over a wide range of magnetic fields {B = 18 
to 34 T), two distinct magnetopolaron branches separated by as much as 11 meV 
(~ 0.4ct;Lo) at resonance (Fig. [12]). The theoretical curves show the results of calculations 
for coupling to the LO phonons in bulk (3D), sheet (2D) and after correction for the 
quasi-2D character of the system using a = 0.29. The agreement between theory and 
experiment is reasonable for the 3D case, but better for the quasi-2D system, if the 
finite spatial extent of the 2D electron gas in the symmetric planar layer is taken into 
account. 

Cyclotron-resonance measurements performed on semiconductor quantum wells 
with high electron density [931 Elj reveal anticrossing near the TO-phonon frequency 
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Figure 12. The cyclotron resonance position plotted as a function of magnetic field 
for InSc. (Reprinted witlr permission from Rcf. [89]. ©1992 by the American Physical 
Society.) 



rather than near the LO-phonon frequency. In Ref. [9l], this effect is interpreted by 
invoking mixing between magnetoplasmons and phonons and in terms of a resonant 
couphng of the electrons with the mixed magneto-plasmon-phonon modes. 

5.2. Cyclotron resonance in a quasi-2D many-polaron system and the role of screening 

In Refs. [95l [96] , the CR spectra for a polaron gas in a GaAs/AlAs quantum well 
are theoretically investigated, taking into account (i) the electron-electron interaction 
and the screening of the electron-phonon interaction, (ii) the magnetoplasmon-phonon 
mixing, (iii) the electron-phonon interaction with all the phonon modes specific for the 
quantum well under investigation. As a result of this mixing, different magnetoplasmon- 
phonon modes appear in the quantum well, which give contributions to the CR spectra. 

It is clearly seen from Fig. [131 that for a high-density polaron gas, anticrossing of 
the CR spectra occurs near the GaAs TO-phonon frequency uti rather than near the 
GaAs LO-phonon frequency uli for both the experimental and the calculated spectra. 
This effect is in contrast with the cyclotron resonance of a low-density polaron gas 
in a quantum well, where anticrossing occurs near the LO-phonon frequency. The 
appearance of the anticrossing frequency close to uti instead of ljli is due to the 
screening of the electron-phonon interaction by the plasma vibrations. A similar effect 
appears also for magnetophonon resonance: as shown in Ref. [9^, the magnetoplasmon- 
phonon mixing leads to a shift of the resonant frequency of the magnetophonon 
resonance in quantum wells from Q ^ u^i to Q ^ ujti- 
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5(T) 

Figure 13. Density map of tlie magnetoabsorption spectra for a 10-nm GaAs/AlAs 
quantum well as calculated in Ref. [25]. Symbols indicate peak positions of the 
experimental spectra (which are taken from Fig. 3 of Ref. [M]). Dashed lines show 
LO- and TO-phonon energies in GaAs. (Reprinted with permission after Ref. [55] . 
©2003 by the American Physical Society.) 

6. Polarons in quasi-OD structures 

6.1. Many-polaron systems in quantum dots 

The ground-state energy and the optical conductivity spectra for a system with a finite 
number of interacting arbitrary-coupling large polarons in a spherical quantum dot 
were calculated using the path- integral formalism for identical particles [981 HOOj - 
A parabolic confinement potential characterized by the confinement energy MIq and 
with a background charge is considered. Using a generalization of the Jensen-Feynman 
variational principle [Ml HOP] , the ground-state energy of a confined A^-polaron system 
is analyzed as a function of N and of the electron-phonon coupling strength a |10HI102] . 

The total spin S is analyzed as a function of the number of electrons in the quantum 
dot for different values of the confinement energy, of the coupling constant and of the 
ratio 7] of the high-frequency and the static dielectric constants. 

Confined few-electron systems, without electron-phonon interaction, can exist in 
one of two phases: a spin-polarized state and a state obeying Hund's rule, depending on 
the confinement frequency (see, e. g., Ref. |103j ). For interacting few-polaron systems 
with « > 3, besides the above two phases, a third phase may occur: the state with 
minimal spin. 

In Ref. |102] the memory-function approach has been extended to a system of 
arbitrary-coupling interacting polarons, confined to a parabolic confinement potential. 
The applicability of the parabolic approximation is confirmed by the fact that the self- 
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induced polaronic potential, created by the polarization cloud around an electron, is 
rather well described by a parabolic potential whose parameters are determined by a 
variational method. For weak coupling, our variational method is at least of the same 
accuracy as the perturbation theory, which results from our approach at a special choice 
of the variational parameters. For strong coupling, the interplay of the electron-phonon 
interaction and the Coulomb correlations within a confinement potential can lead to 
the clustering of polarons in multi-polaron systems. Our choice of the model variational 
system is justified because of this trend, which occurs in a many-polaron system with 
arbitrary A^, for a finite confinement strength. 
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Figure 14. Optical conductivity spectra of = 20 interacting polarons in CdSe 
quantum dots with a ~ 0.46, rj ~ 0.656 for different confinement energies close to 
the transition from a spin-polarized ground state to a ground state obeying Hund's 
rule. Inset: the first frequency moment (f2) of the optical conductivity as a function 
of the confinement energy. (Reprinted with permission from Ref. |102j . ©2004 by the 
American Physical Society.) 



The shell filling schemes of an A^-polaron system in a quantum dot can manifest 
themselves in the optical conductivity. In Fig. [HI optical conductivity spectra for 
= 20 polarons are presented for a quantum dot with the parameters of CdSe: 
a = 0.46, rj = 0.656 [104j and with different values of the confinement energy ^JIqEI In 

* For the numerical calculations, we use effective atomic units, where h, the electron band mass 
and e/y/Scxi have the numerical value of 1. This means that the unit of length is the effective Bohr 
radius a*g = h^£oo/ (w&e'^), while the unit of energy is the effective Hartree H* = Jn^e^/ (fi^e^). 
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this case, the spin-polarized ground state changes to the ground state satisfying Hund's 
rule with increasing WLq in the interval 0.0421iJ* < MIq < 0M22H*. 

In the inset to Fig. [HI the first frequency moment of the optical conductivity 
/~ QRea (Q) dQ 



r Rea (n) dn 



(16) 



as a function of HQo shows a discontinuity, at the value of the confinement energy 
corresponding to the change of the shell filling schemes from the spin-polarized ground 
state to the ground state obeying Hund's rule. This discontinuity could be observable 
in optical measurements. 

The shell structure for a system of interacting polarons in a quantum dot is clearly 
revealed when analyzing the addition energy and the first frequency moment of the 
optical conductivity in parallel. In Fig [151 we show both the function 

e(iv) = (fi)U^i-2(fi)u+ (17) 

and the addition energy 

A (A^) =E^{N + 1)- 2E° (N) +E'^{N-1). (18) 
for interacting polarons in a 3D CdSe quantum dot. 
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Figure 15. The function 8 {N) and the addition energy A (A^) for systems of 
interacting polarons in CdSe quantum dots with a = 0.46, 77 = 0.656 for Mlf) = Q.IH* . 
(Reprinted with permission from Ref. |102| . ©2004 by the American Physical Society.) 



As seen from Fig [TJl distinct peaks appear in G (A^) and A (A^) at the "magic 
numbers" corresponding to closed-shell configurations at = 8, 20 and to half-filled- 
shell configurations at A^ = 5, 14. We see that each of the peaks of 6 (A^) corresponds to 
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a peak of the addition energy. The peak positions for the addition energy of interacting 
polarons in a 2D parabohc quantum dot |102] agree well with the experimental results 
for the addition energies of cylindrical GaAs quantum dots |106] . 

6.2. Non-adiahaticity of polaronic excitons in semiconductor quantum dots. 
Photoluminescence and Raman scattering of polarons in quantum dots 

Interest in the optical properties of quantum dots has been continuously growing i. a. 
because of the prospects of these structures for optoelectronic applications. Recent 
measurements [1071 11081 11091 [TTO] of photoluminescence of self-assembled InAs/GaAs 
quantum dots, reveal unexpectedly high probabilities of the phonon-assisted transitions. 
Attempts to interpret some of these experiments on the basis of the adiabatic theory 
meet with difficulties. For spherical CdSe quantum dots, e. g., the Huang- Rhys factor 
S calculated within the adiabatic approximation takes values which are significantly 
(by one or two orders of magnitude) smaller than the so-called "experimental Huang- 
Rhys factor". This "experimental Huang- Rhys factor" is determined from the ratio of 
the measured intensities of the phonon satellites. In the framework of the adiabatic 
approach, various mechanisms, which cause the separation of the electron and hole 
charges in space |108l lllll I112j , are commonly considered as a possible origin for the 
unexpectedly large Huang-Rhys factor. 

It has been shown in Ref. |105j . that non-adiabaticity of exciton-phonon systems in 
some quantum dots drastically enhances the efficiency of the exciton-phonon interaction, 
especially when the exciton levels are separated with energies close to the phonon 
energies. Also "intrinsic" excitonic degeneracy can lead to enhanced efficiency of the 
exciton-phonon interaction. The effects of non-adiabaticity are important to interpret 
the surprisingly high intensities of the phonon 'sidebands' observed in the optical 
absorption, the photoluminescence and the Raman spectra of some quantum dots. 
Major deviations of the observed oscillator strengths of some pho non-peak sidebands 
from the Franck-Condon progression, which is prescribed by the commonly used 
adiabatic approximation, find a natural explanation within our non-adiabatic approach, 
introduced in [105], see [ml [IH fTT5lfTT6] . 

In Ref. |105j . a method was proposed to calculate the optical absorption spectrum 
for a spherical quantum dot taking into account the non-adiabaticity of the exciton- 
phonon system. This approach has been further refined in Ref. |117] : for the matrix 
elements of the evolution operator a closed set of equations has been obtained using a 
diagrammatic technique. This set describes the effect of non-adiabaticity both on the 
intensities and on the positions of the absorption peaks. 

While the adiabatic approximation predicts negligibly low intensities of one- and 
two-phonon sidebands, our non-adiabatic theory allows for a quantitative interpretation 
of the observed high intensity of the phonon sidebands in the photoluminescence 
(Fig. [T6l) and Raman (Fig. [T71) spectra of some quantum dots, in agreement with 
experiment. 
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Figure 16. Photoluminescence spectra of colloidal spherical CdSc quantum dots with 
wurtzite structure. The dashed curve represents the experimental data from Ref. [118j . 
The dash-dotted curve displays a result of the adiabatic approximation - a Franck- 
Condon progression with Huang-Rhys factor S — 0.06 as calculated by S. Nomura and 
S. Kobayashi [lllj . The solid curve results from the non- adiabatic theory. (Reprinted 
with permission after Ref. [105j . ©1998 by the American Physical Society.) 
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Figure 17. Resonant Raman scattering spectra of an ensemble of CdSe quantum 
dots with average radius 2 nm at T = 77 K (panel a) and of PbS quantum dots with 
average radius 1.5 nm at T = 4.2 K (panel h). The dash-dot-dot curves show the 
luminescence background. The dash-dot curves in panel b indicate contributions of 
phonon modes (classified in analogy with electron states in a hydrogen atom) into 
the Raman spectrum. The dashed curves represent experimental Raman spectra. 
(Reprinted with permission after Ref. |114| . ©2002 by the American Physical Society.) 
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From the comparison of the spectra obtained in the adiabatic approximation 
with those resulting from the non-adiabatic approach, the following effects of non- 
adiabaticity are revealed. First, the polaron shift of the zero-phonon lines with respect 
to the bare-exciton levels is larger in the non-adiabatic approach than in the adiabatic 
approximation. Second, there is a strong increase of the intensities of the phonon 
satellites as compared to those derived within the adiabatic approximation. Third, 
in the optical absorption spectra found within the non-adiabatic approach, phonon 
satellites appear, that are related to non-active bare exciton states. Fourth, the optical- 
absorption spectra demonstrate the crucial role of non-adiabatic mixing of different 
exciton and phonon states in quantum dots. This results in a rich structure of the 
absorption spectrum of the exciton-phonon system |119[ I113[ I115j . Similar conclusions 
about the pronounced influence of the exciton-phonon interaction on the optical spectra 
of quantum dots have been subsequently formulated in Ref. [120] in terms of a strong 
coupling regime for excitons and LO phonons. Such a strong coupling regime is a 
particular case of the non-adiabatic mixing, studied by us [10511119] . related to a (quasi-) 
resonance which arises when the spacing between exciton levels is close to the LO phonon 
energy. The large enhancement of the two-phonon sidebands in the luminescence spectra 
as compared to those predicted by the Huang- Rhys formula, which was explained in Ref. 
|105[ 1119] by non-adiabaticity of the exciton-phonon system in certain quantum dots, 
has been reformulated in Ref. |120j in terms of the Frohlich coupling between product 
states with different electron and/or hole states 

Due to non-adiabaticity, multiple absorption peaks appear in spectral ranges 
characteristic for phonon satellites. From the states, which correspond to these 
peaks, the system can rapidly relax to the lowest emitting state. Therefore, in the 
photoluminescence excitation (PLE) spectra of specific quantum dots, pronounced peaks 
can be expected in spectral ranges characteristic for phonon satellites. New experimental 
evidence of the enhanced phonon-assisted absorption due to effects of non-adiabaticity 
has been provided by PLE measurements on single self- assembled InAs/GaAs quantum 
dots [121] and InGaAs/GaAs quantum dots |122j . 
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